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1. INTRODUCTION 

Satellites are exposed to many disturbances in space [1] which can adversely affect their attitude. 
The reliability of the microsatellite data acquisition is highly dependent on the quality of the SAC system 
because the attitude control system is crucial for accurate satellite operation [2]. Therefore, the needs for 
an accurate design of attitude control subsystem which guarantees the quality and reliability of the data. 
Several methods exist in the literature for achieving SAC. [3] and [4] proposed direct and modified adaptive 
controllers, respectively to mitigate different scenarios in spacecraft attitude control. Fixed-time [5], 
finite-time [6], adaptive gain based integral second order [7], fast terminal [8] and adaptive fast terminal [9] 
sliding mode techniques have been employed in rigid spacecraft attitude controller design. 

Several applications of PID controller for satellite attitude control can be found in the literature. PID 
attitude controller for a rigid body satellite designed using a combination of Kharitonov theorem 
and vectored particle swarm optimization method was presented in [10]. Spacecraft control torque 
and velocity constraints dependent sliding mode based robust PID attitude controller were presented 
by Y. Li et al. [11]. Latest research by [12-16] employed an observer-based PID controller to control 
the attitude of a rigid satellite with ON/OFF thruster actuator which showed improvement in attitude control 
performance. The optimal PID controller for SAC, designed by [17] gave acceptable results for both 
the settling time and percentage overshoot, even though the bounds of controller parameters were arbitrarily 
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chosen [18] designed both PID and Fuzzy logic-based satellite attitude controllers in which the dynamics 
model of BILSAT-1 were utilized. Compared to the PID controller, the Fuzzy logic counterpart performed 
better in terms of system response speed and robustness. More so, [19] designed a PID controller to optimally 
control the satellite yaw attitude; and when tested via simulation, the settling time and percentage overshoot 
obtained were satisfactory as their values were within the design specifications. 

PID controllers are popularly used in the control of industrial processes due to their structural 
and functional simplicity which in turn make their tuning relatively simple, most especially by a trial 
and error method [20]. It is also found in the majority of control systems due to its reliable performance and 
fine control capabilities [21]. At least 95% of control loops in process industries are PID in nature [22], 
which can be attributed partly to their robust performance in a wide range of operating conditions [23]. 
The performance of PID controller can be greatly improved by using an ITAE robust controller design which 
optimizes the controller gains, and subsequently the system performance parameters. The tool was adopted in 
this study as the system optimum performance parameters can be easily extracted as the system parameters 
and controller gains are varied [23]. Also, ITAE was identified as one of the best methods for robust control 
system designs [24]. The optimum coefficients for the ITAE performance criterion for a step input that can 
be used to find the best controller parameters have been determined for the standardized closed-loop transfer 
function in [23, 24]. Compared the results obtained by using Ziegler-Nichols and ITAE tuning methods for 
both first-order and second-order systems with dead time [25]. The results showed that ITAE tuning method 
was better than the Ziegler-Nichols tuning method in both scenarios. In most cases, the closed-loop response 
of the controlled system may not meet the desired specifications, under this condition a pre-filter transfer 
function must be designed and connected in series with the designed closed-loop system transfer 
function [24, 26]. It should be noted that all the parameters of the designed pre-filter transfer function are 
normally obtained directly from the dc gain and zeros of the designed closed-loop system transfer 
function [26]. 


2. DESIGN SPECIFICATIONS 

The system design specifications considering unit step input are shown in Table 1. In this study, 
the system PO is a function of system damping ratio, while Ts is the time for which 
the system response remains within 2% of the final value as given in (1) and (2), respectively. 


Table 1. Design specifications 











Parameters Values 
Percentage Overshoot (PO) <5% 
Settling time, T, (for 2% criterion) <2s 
Steady-state error 0 
P, = 100e ("V7") (1) 
4 
Tee (2) 


Where Cand w, are the damping ratio and natural frequency in rad/sec for the closed-loop system. From (1): 


X2 
= Jax p 


where X = —log0.01P0. Also, from (2), 


ona (4) 


in order to achieve the system design specifications, the threshold values for ¢ and w, were determined 
using (3) and (4) [27]. The required conditions are ¢ > 0.7797 (~0.8) and w, = 6 rad/s. 
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3. SYSTEM MODEL 

The yaw-axis attitude model for the satellite used in this study was derived based on a single input 
single output and linear time-invariant assumption which then yielded the transfer function that was used for 
the controller design. 


3.1. Determination of open-loop transfer function and PID control algorithm selection 
The yaw-axis ACS is represented by the block diagram of Figure 1 while its sub-systems’ transfer 


functions are as expressed in (5)-(8). 








8,(s) a V.(s) G(s) V(s) Gamp(8) v,(s) aa ws 




















8,;(s) 
Figure 1. Block Diagram of microsatellite yaw-axis ACS 
= 6,(s) _ 1 

Gsat(S) = LG) ~ 0.882 (5) 

_ T,(s) _ 78.38 
Gld) = 7 (5) = 52 4 1815.45 + 24466 © 

O W(s)_ 240 

Gamel) = We) OAs +1 o) 

_ V(s) _ Kps? + Kps +K, 
G6) =p = (8) 


Where all the parameters are in Laplace Transform (s-domain); 6, is the actual yaw attitude in degrees; 
T, is the actuator torque in Nm; V, is the motor armature voltage in volts; V is the controller output voltage 
in volts; V, is error voltage in volts; G,,, is the satellite structure dynamic transfer function; G,(s) 
is the actuator dynamic transfer function; Gamp(S) is the amplifier transfer function; G,(s) is the controller 
transfer function; K, is the proportional gain; K; is the integral gain; and Ką is the derivative gain. 
The system open-loop transfer function shown in (9) was determined by cascading the transfer function of (5) 
to (8). The closed-loop transfer function obtained from (9) is as expressed in (10). 


B(s) _ 234900 * Kas? + 234900 * K,s + 234900 « K; 7 
V,(s) s5 + 1825s* + 426255? + 24462582 (9) 





6,(s) 234900 * Kas? + 234900 * K,s + 234900 + K; 
0,(s) s5 +1825 *s* + 42625 *s3 + (244625 + 234900 * K4) * s? + 234900 * K,s + 234900 « K; 








(10) 


3.2. Reduction of system transfer function for controller design 

An analysis of the uncontrolled system response to a specified test signal is necessary to determine 
the degree of the control effort required if any. To achieve this, the uncontrolled yaw-axis ACS open 
and closed-loop transfer functions, (11) and (12), were obtained from (9) and (10) respectively, assuming 
a unity controller transfer function. (12) was subjected to unit step function input and the resulting system 
step response is as shown in Figure 2. As shown in the Figure, the response will reach the set point at infinite 
time, therefore there is no overshoot, and settling time was 3.49 seconds. Considering this step response there 
is the need to incorporate a controller into the system so that the system design specifications shown in 
Table 1 can be achieved. 








By(s) _ 234900 di 
V(s) s4 +1825s3 + 42625s? + 244625s ) 
Oy(s) _ 234900 2) 
6,(s)  s*+1825s3 + 42625s2 + 244625s + 234900 
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The PID control algorithm was considered due to its numerous advantages, such as simple design 
methodology and effectiveness in system control. It should be noted that PID control algorithm (either in full 
PID, PI or PD) performs better when the order of the system to be controlled is first, second or higher-order 
that can be approximated by either first or second-order dynamic [23]. The system to be controlled was 
reduced from 4" to a 2" order using Balred reduction model, and the transfer function of the reduced-order 
version is as shown in (13). Both the 4" order and 2™ order versions exhibit good conformity, as evidenced 
in the plots of unit step responses of (12) and (13) shown in Figure 3. 
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Figure 3. Plots of closed-loop unit step responses of the 4" and 2"4 order of the uncontrolled system 


The step response parameters of the two systems shown in Table 2 were obtained from Figure 3. 
Considering the closeness of the values, the higher-order one can be approximated by its reduced order 
version, which is subsequently used for PID controller design to achieve better performance. 


Table 2. Step response parameters for uncontrolled higher and reduced-order closed-loop systems 





Time-domain parameters 


Higher-order system _ Reduced-order system 





Rise time (sec.) 
Settling time (sec.) 
Peak amplitude 
Peak time (sec.) 


% Overshoot 


Steady state final value 


1.890 1.800 
3.490 3.290 
0.998 0.997 
5.230 5.580 
0.000 0.000 
1.000 0.998 
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Grou (s) 


Gre(S) = IF GaG) (14) 


Substituting (13) into (14) yields the open-loop transfer function of the reduced-order system Go: (S), 
as follows: 


—0.16728s + 6.081 


_ _—0.16728s + 6.081 _ 15 

Gro(S) = 34616938 + 0.0150 = 
As +B 

Gro(S) = s2 + Cs +D S 


where A = -0.1673, B = 6.081, C = 6.1693 and D = 0.0150. 


4. CONTROLLER DESIGN 
4.1. Determination of PID controller gains 

The full PID controller transfer function is as expressed in (8), while the open-loop PID controlled 
transfer function for the reduced-order system was obtained by multiplying (8) and (16) to yield (17). 





AKas°? + (BK4 + AKy)s? + (BK, + AK;)s + BK; 
s3 + Cs? + Ds 


Grcol (s) = ( (17) 


Using (17) the closed-loop transfer function for the reduced-order PID controlled system, G,¢,)(s) is 
given by (18). 


(ye) (AK,s? + (BK; + AK,)s? + (BK, + AK;)s + BK;) 

















Greei(s) = , . (BKa+AKp +0), , (BKp t+ AK +D BK, (18) 
s +( T+AK, ) +( T+AKyz Js +r ie 
The characteristic of (18) is (19). 
: (Gece A A , BK _ (19) 
1+AK, 1+AK, 1+AK, 


Since the characteristic equation of the transfer function is of 3" order, the value of controller gains 
using ITAE design method can be determined by obtaining a standardised form of 3" order characteristic 
equation. These can be found in the table of optimum coefficients of system closed-loop transfer function 
based on the ITAE criterion for a step input [23]. The obtained standardised characteristic equation from 
the table is given in (20). 


s? + 1.75wps? + 2.150,2s + w,? = 0 (20) 


Where w,= 6 rad/s from (4). Equating the coefficients of s?,s* and s in (19) and (20), and then rearranging, 
we can obtain 


—0.1673K, + 7.8377Kq = 4.3307 (21) 
6.081K, — 0.1673K; + 13Kq = 77.385 (22) 
K; + 5.9426Kq = 35.5205 (23) 


The values of Kp, K; and Kg are 11.8492, 30.7337, and 0.8055 respectively, obtained by 
solving (21), (22) and (23) simultaneously. The transfer function G,,.; of the closed-loop PID-controlled 
reduced-order system was obtained by substituting the values of A, B, C, D, Kp, Kj and Kg in (18). 


—0.1557s3 + 3.37s? + 77.335 + 216 


24 
s3 + 10.5s 2 + 77.35s + 216 ( ) 


Greet = 
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Also, the transfer function Ghee) of the higher-order system of the PID was obtained by substituting 
the values of Kp, K; and Kg in (10). 


189200(s? + 14.71s + 38.16) 


Gneel = 25 
heel “"'s5 + 1825s4 + 42625s3 + 433800s? + 2783000s + 7219000 a 





4.2. Determination of PD controller gains 
The uncontrolled system reduced-order open-loop transfer function was already given in (16). 
The PD controller transfer function was obtained from (8) by assuming K; = 0 and is as expressed in (26). 


Ge = (Kas + Kp) (26) 
The open and closed-loop PD controlled reduced-order system transfer function were derived based on (16) 


and (26), as shown in (27) and (28) respectively. 
_ (AKas? + (BKa + AK, )s + BKy) 

















Greot s2 +Cs + D ep 
1 
E Gs AK) (AKys? + (BKa + AKy)s + BK) 
reel > , (BKat AKyp +C. BK, +D (28) 
s ( I FAK, ) I+AK; 
The characteristic (28) is as shown in (29). 
F + + AK, + e BK, +D sb 
z 1+AK 1+AKı 2) 


The characteristic of (29) is of order 2, therefore to determine the PD controller gains values by 
applying ITAE method the corresponding standardised second-order characteristic equation was obtained 
from the optimum coefficients’ table. The obtained characteristic equation is: 

s? + 1.40,8 + Wy” (30) 
where W,= 6 rad/s. Equating the coefficients of s* and s in (29) to (30) and rearranging, we have 
—0.1673Ky + 7.4863Kq = 2.2307 (31) 


6.081K, + —6.0228Kq = 35.9850 (32) 


the values of K, and Kg were obtained as 5.5008 and 0.4209, respectively, by solving (31) and (32) 
simultaneously. Substituting the values of A, B, C, D, Kp and Kg in (28) gives (33). 


_ 0.07575(—s? + 23.27s + 475) 
ie (s2 + 84s + 36) 





(33) 


Substituting the values of Kp, Ka with K;=0 into (10) gives the closed-loop transfer function of PD controlled 
higher-order system of (34) 


98870(s + 13.07) 


a o a 34 
s* + 1825s3 + 42625s2 + 343500s + 1292000 ( ) 


Gheel = 


5. DESIGN OF PRE-FILTER TRANSFER FUNCTION 

As can be seen from (24), (25), (33) and (34) both the reduced and higher-order PID and PD 
controlled closed-loop systems have zeros, which will negatively affect their response [23, 24, 26]. 
To ameliorate this effect, there is the need to incorporate pre-filter transfer function in series with the system 
closed-loop transfer function as shown in Figure 4. This is necessary to ensure that the resulting closed-loop 
transfer function will contain no zero and still maintain unity dc gain which will lead to the overall 
improvement of the system response [23]. Two pre-filter transfer functions were designed, the first 
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was based on the reduced-order system and the second was based on the higher-order system for each of 
the controller design. 


8,(s) 


— 





Figure 4. System closed-loop transfer function with pre-filter 


5.1. Determination of pre-filter transfer function using PID controlled reduced-order system 
Factorizing the numerator of (24) yielded (35). 


0.1557(s — 36..3250)(s + 11.3653) (s + 3.3596) 


mam s3 + 10.582 + 77.35s + 216 (353) 

It can be seen from (35) that the introduction of PID controller transfer function into the system has added 2 
more zeros to this system at s = —11.3653 and s = —3.3596 in addition to the already existing zero for 
the uncontrolled system which is at the point s = 36.3250. The expected closed-loop reduced-order PID 
controlled system transfer function when the pre-filter transfer function is introduced is as given in (36). 


216 


G = 36 
reclf ~ 53 + 10.582 + 77.358 + 216 (36) 





But, 
Grecif = Greet Grf (37) 


therefore, the lower order-based pre-filter transfer function (G,-) was realized using (37) and is as 
expressed as 


a 1387 
"f (—s3 + 21.6s? + 496.7s + 1387) 





(38) 


As can be seen from (38), one of the poles of the reduced-order pre-filter transfer function is on 
the right-hand side of s-plane, therefore combining it in series with the controlled system closed-loop transfer 
function will cause instability. As a result, the response of the closed-loop controlled system with pre-filter 
will be different from the expected response of (36). In order to improve the system response, obtaining 
the pre-filter transfer function from the closed-loop PID controlled transfer function for the original 
higher-order system was considered. 


5.2. Determination of pre-filter transfer function using higher-order PID controlled system 

The PID controlled closed-loop transfer function for the original higher-order system is as 
shown in (25). Factorizing the numerator of (25) yielded (39), which shows that incorporation of PID 
controller into the system has increased the number of system zero from 0 to 2 and these zeros are at 
s = —11.3470 and s = —3.3630. 


189200(s + 11.3470) (s + 3.3630) 


5 4 3 2 (39) 
5° + 1825s* + 42625s° + 433800s* + 2783000s + 7219000 


Gheel = 


Now the transfer function for the controlled closed-loop system with a pre-filter transfer function 
(Ghear) is as expressed in (40). 


Gneat = Gre Geel (40) 
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Where: Ghr is the pre-filter transfer function obtained from the original higher-order control system. 
The expected controlled system transfer function with pre-filter is as expressed in (41). 


Be pr 7219000 (41) 
heclf — S5 + 1825s4 + 42625s? + 433800s2 + 2.783000s + 7219000 
Therefore, from (39), (40), and (41), 
38.16 
Gn (42) 


~ (s2 + 14.71s + 38.16) 


Since Gp, is stable and can be used to improve the system response, therefore the transfer function 
for PID controlled higher-order system with higher-order pre-filter using (40) can be expressed as 





38.16 ) ( 189200(s? + 14.71s + 38.16) ) 


Great = ( 
heclf — \s2 + 14.71s + 38.16/ \s5 + 1825s* + 42625s? + 433800s2 + 2783000s + 7219000 


(43) 


5.3. Determination of pre-filter transfer function from reduced-order PD controlled system 
The PD controlled reduced-order closed-loop system transfer function is as shown in (33), which 
may be rewritten as, 


0.07575(s — 36.3407)(s + 13.0707) 
recl 7 ae a Aa: E SAA = (44) 
(s* + 8.4s + 36) 
This shows that the introduction of PD controller into the system has increased the number of reduced-order 
system closed-loop zeros from 1, for the uncontrolled system located at s = 36.3522 to 2 for the controlled 
system at s = 36.3522 ands = —13.0707. The expected closed-loop reduced-order PD controlled system 
with pre-filtered transfer function is as given in (45). 


36 


es ee 4 
Greet = 24 B45 + 36 oe) 
But from (37), 
G 
Grp = (46) 
reclf 
therefore, 
475.2475 
Grf = (47) 





(—s? + 23.27s + 475) 


The closed-loop reduced-order PD controlled system with this pre-filter transfer function will not 
give the expected pre-filtered closed-loop response because of the introduction of the right-hand side pole 
into the system by the pre-filter transfer function. Since the pre-filter transfer function obtained from 
the closed-loop reduced-order PD controlled system will destabilize the system, the next option is to obtain 
the pre-filter transfer function from the closed-loop PD controlled transfer function for the original 
higher-order system. 


5.4. Determination of pre-filter transfer function from higher-order PD controlled system 
The PD controlled closed-loop transfer function for the original higher-order system is as 
shown in (34). Factorizing the numerator of (34) gives (48). 


98870(s + 13.07) 


Cee 48 
hccl — <4 + 1825s3 + 42625s2 + 343500s + 1292000 (48) 


It can be seen from (48) that the introduction of PD controller into the system has increased the number of 
reduced-order system closed-loop zero from 0 for the uncontrolled system to 1 for the controlled system at 
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s = —13.0707. The expected PD controlled closed-loop transfer function for the original higher-order 
system with pre-filter transfer function is as shown in (49). 








Be 1292000 (49) 
hcclf — s4 + 182583 + 42625s2 + 343500s + 1292000 
But from (40), 
Ghccl 
G= (50) 
Gheclt 
therefore, 
ne 13.07 = 
hf — (s + 13.07) 61) 


From (51), the transfer function for PD controlled higher-order system with higher-order pre-filter is 
expressed as, 


13.07 98870(s + 13.07) ) 


ee 52 
s+ 307) (s + 1825s? + 42625s2 + 343500s + 1292000 ( ) 


Gheelf = ( 


6. SIMULATION 

The controller design carried out in this paper has resulted into four types of controlled systems, 
namely PID controlled system without pre-filter, PID controlled system with pre-filter, PD controlled system 
without pre-filter and PD controlled system with pre-filter. It should be noted that the system of interest 
is the higher-order system which is the original system. Therefore, to select the most appropriate out of these 
four designed micro-SAC systems, there is the need to investigate their performance via simulation. 
Therefore, the system unit step response was plotted for these systems and the time domain performance 
parameters were extracted and used as the basis for selecting the most appropriate system after comparison. 
The time-domain performance parameters of interest are the rise time (T,), settling time (T,), percentage 
overshoot (PO), time to peak overshoot (Tp) and steady-state error (E,,). To easy the system performance, 
comparison and analysis the system simulation was carried out using the uncontrolled and controlled systems 
closed-loop transfer functions as follows. Using (12), (25), (34), (43) and (52), the closed-loop step responses 
for the uncontrolled higher-order system; PID controlled higher-order system without pre-filter; PD 
controlled higher-order system without pre-filter; PID controlled higher-order system with pre-filter and PD 
controlled higher-order system with pre-filter were plotted as shown in Figure 5. 











Uncontrolled higher order system 
PID controlled higher order system without pre-filter 
PID controlled higher order system with pre-filter 
— — — PD controlled higher order system without pre-filter 
— — — PD controlled higher order system with pre-filter 

1 f 1 








Satellite yaw axis attitude (degrees) 











1 2 3 4 5 
Time (seconds) 


Figure 5. Closed-loop step response for uncontrolled and controlled higher order system 
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7. RESULTS AND DISCUSSION 

The step response parameters in Table 3 were obtained from the uncontrolled and controlled 
systems closed-loop step response of Figure 5. As can be seen from Figure 5 and Table 3, the PID-controlled 
system without pre-filter has a better performance compared to the uncontrolled system when the rise time, 
settling time and time to peak overshoot were considered, but the performance becomes worse in terms 
of percentage overshoot. Introduction of pre-filter into the PID controlled system reduces the settling time 
and percentage overshoot but worsen both the rise time and time-to-peak overshoot by increasing the 
numerical value of these parameters. Introduction of pre-filter into the PD controlled system worsen the 
system performance by increasing the numerical value of the rise time, settling time and time to peak 
overshoot. But there is a reduction in the value of the percentage overshoot encountered by this system when 
compared to PD controlled system without pre-filter. Both the PD controlled systems without and with pre- 
filter have better performance compared to the uncontrolled system considering the rise time, settling time 
and time to peak overshoot. But the percentage overshoot for the uncontrolled system is less than that for any 
of the two PD controlled system. 


Table 3. Step response performance parameters for yaw-axis attitude closed-loop controlled system 








System condition T, (s) T; (s) PO Tp (s) Ez 
Uncontrolled system 1.89 3.49 0 5.23 0 
PID controlled system without pre-filter 0.136 1.31 48.1 0.37 0 
PID controlled system with pre-filter 0.365 1.17 1.07 0.761 0 
PD controlled system without pre-filter 0.288 0.814 4.73 0.593 0 
PD controlled system with pre-filter 0.332 0.887 3.66 0.706 0 





From the results, there are three major time-domain performance specifications to be satisfied; it can 
be seen in Table 3 that the uncontrolled system satisfies the PO and E;, conditions, PID controlled system 
without pre-filter satisfied the T, and E,, conditions, PID controlled system with pre-filter satisfied 
the T,, PO and Ess conditions, the PD controlled system either with or without pre-filter satisfies the Ts, 
PO and E,s conditions. Finally, three out of the four controlled systems designed satisfy the design 
specifications of T, <2 s, PO < 5%, and Ess = 0. It is noted that the three are the PID controlled system with 
pre-filter, PD controlled system without pre-filter and PD controlled system with pre-filter. Based on this 
analysis any of the three systems can be used but PD controlled system without pre-filter is the most 
preferable due to its structural simplicity which will enhance its digital implementation. 


8. CONCLUSION 

PID and PD controller-based microsatellite attitude controllers were designed. Pre-filter transfer 
functions were designed and incorporated into the designed systems to improve their performance. 
The designed systems were tested via simulation and the results revealed that three out of the four systems 
satisfy the design specifications and the PD controlled system without pre-filter transfer function is the most 
preferable due to its structural simplicity which will enhance its digital implementation. 
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